Trapping in quantum chains 
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Abstract 

A quantum breather on a translationally invariant one-dimensional 
anharmonic lattice is an extended Bloch state with two or more particles 
in a strongly correlated state. We discuss several effects that break the 
lattice symmetry and lead to spatial localization of the breather. 

1 Introduction 

The localization of energy by nonlinearity in classical lattices has been much 
studied in the last 20 years. The corresponding localized states, known as 
intrinsic localized modes or discrete breathers, have been the subject of intense 
theoretical and experimental investigation pQ. Corresponding results on the 
quantum equivalent of discrete classical breathers are less numerous, c.f. |3 |3] 
for some theoretical results and 0] for some experimental work. Studies of 
quantum modes on small lattices may be relevant to studies of quantum dots 
and quantum computing (c.f. 

In this paper, we present some results related to quantum lattice problems, 
in particular in one dimensional lattices with a small number of quanta. We 
study a periodic lattice with / sites containing bosons, described by the quan- 
tum discrete nonlinear Schrodinger equation (QDNLS), a quantum version of 
the discrete nonlinear Schrodinger equation (also know as the Boson Hubbard 
model). The DNLS equation describes a particularly simple model for a lattice 
of coupled anharmonic oscillators, and it has been used to describe the dynam- 
ics of a great variety of systems [H]. The corresponding quantum Hamiltonian 
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is given by 



H = - 



t 1 



~ fj btbtb j b j +e j bt(b j -i+b j+1 ), 



(1) 



where and 6j are standard bosonic operators, Jj/ej is the ratio of anhar- 
monicity to nearest neighbor hopping energy, and the chain is subject to peri- 
odic boundary conditions with period /. Initially we consider the case where 
the chain is translationally invariant, i.e. 7j = 7 and tj — e are independent of 
j. In general we take e = 1. 

The Hamiltonian has an important conserved quantity, the number N = 



J2j=i bjbj, which enables the total Hamiltonian to be block-diagonalised and 
greatly simplifies the analysis. In this letter we restrict ourselves to the simplest 
nontrivial case N — 2, though many of the results are valid for larger values 
on N. The N — 2 bound states corresponds to bound two-vibron states, as 
observed experimentally in several systems [7]. 

In QDNLS case, we use a number state basis, \ip n ) = [ni, 712, n/], where rii 
represents the number of quanta at site i (N = Y] m). A general wave function 

is |*„) = ^2 n Cn\lp n }- 

In homogeneous quantum lattices with periodic boundary conditions, it is 
possible to block-diagonalize the Hamiltonian operator using eigenfunctions of 
the translation operator with fixed value of the momentum k [2]- 



Figure 1: Eigenvalues E(k) for QDNLS model with periodic boundary condi- 
tions. N = 2, / = 19, 7 = 4. 



As shown in Fig. ^ if the anharmonicity parameter is high enough, there 
exists an isolated eigenvalue for each k which corresponds to a localized eigen- 
function. By this we mean there is a high probability of finding the two quanta 
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on the same site, but due to the translational invariance of the system, an equal 
probability of finding these two quanta at any site of the system. In these cases, 
some analytical expressions can be obtained in some asymptotic limits (solutions 
of this problem go back to the 30's, for recent discussions see [51 121 OH )■ 

In particular, working at k — for simplicity, the ground state unnormalized 
eigenfunction is 

|*) = [20 ... 0] + [020 . . . 0] + • • • + [0 . . . 02] + Oi-f- 1 ), 

i.e. on a lattice of length /, the unnormalized coefficients c,; of the first / terms 
are equal to unity and the rest are 0(7 _1 ). In the other extreme of complete 
spatial localization, one of these Cj would be unity and the rest zero. In this 
letter we consider how these components change as the translational invariance 
of the lattice is broken in various ways. 

One simple way that translational invariance can be broken is by considering 
a finite chain with no-flux boundary conditions. The Hamiltonian operator now 
cannot be block-diagonalizcd using eigenvectors of the translation operator. In 
this case, the computational effort increases, but it is still possible to calculate 
the all the eigenvalues and eigenvectors of the Hamiltonian operator, if / and 
N are small enough, by using algebraic manipulation methods and numerical 
eigenvalue solvers. 

The existence of local inhomogeneities or impurities in a system can affect the 
nonlinear localized modes considerably. Also, in systems with both nonlinearity 
and impurities, it is important to understand the interplay between these two 
sources of localization. For these cases we break translational invariance by 
making one or more of the jj or the €j depend on j. This may occur because of 
localized impurities, or because the chain geometry becomes non- uniform. Two 
examples of non-uniform geometries are shown in Fig. [5] 




Figure 2: Two non-uniform chain geometries 

Fig-EK shows a circular chain twisted into a figure-of-eight, so that two sites 
on the chain, which arc distant measured along the length of the chain, become 
spatially close. This toy model for a globular protein was studied in ,9^, where 
it was shown that moving breathers described by the classical DNLS equation 
could become trapped at the cross-over point. In the quantum case such a 
geometry can be modeled by adding a term such as 

(blb m + blb e ), (2) 
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to the Hamiltonian, where I and m are the two sites brought close together by 
the twist, and ai yTn is the separation distance relative to the unit length of the 
unperturbed chain. This can be considered as a special case of a chain with 
long-range coupling. See also ^U] f° r a more realistic protein simulation, and 
[TT] for other discussions on the effects of chain geometry on moving breathers. 

The bent chain in Fig. \Bp shows another possible geometry which has been 
studied recently in the classical DNLS case. In this case we have an abrupt 
bend which is simulated by adding an additional term as in @ but where 
m = mo — 1,£= too + 1> where mo is the vertex of the bend. By varying the 
values of a^ m = a, all angles between and n can be simulated approximately. 
The influence of this geometry has been analyzed in the DNLS context in |12|. 
and in nonlinear Klein-Gordon systems in ^21- This geometry is of interest in 
nonlinear photonic crystals waveguides and circuits |14j. 

Localization due to random variation of the lattice parameters has long been 
studied in the harmonic model since the pioneering work of Anderson J^l • Our 
interest is to see what new localization effects the anharmonic terms bring to the 
model, and to what extent the anharmonic effects enhance the Anderson-like 
localization effect when this is present in the harmonic model (7, = 0) . See also 
the discussions in [Po| . 

Most of our findings are not specific to the QDNLS model, for example we 
have repeated our calculations using the attractive fermionic Hubbard model 
with two particles of opposite spins. Details will be given elsewhere. Although 
all the models we consider have a conserved number, in general they are not 
quantum integrable. 

We now consider the effects introduced above in more detail. 

2 Localization in an straight chain with impuri- 
ties 

In this version of our model, in order to explore the interplay between the lo- 
calization induced by the nonlinearity and the influence of a impurity in these 
localized states, we introduce a local inhomogeneity in the anharmonic parame- 
ter. To isolate the effect of the impurity of other effects related to the finite size 
of the chain, we retain the periodic boundary conditions. The anharmonicity 
parameter is 7^ = 7j m , and jj = 7 for j ^ i. 

In the homogeneous system (7 im = 7), with 7 large enough, as discussed 
above, the ground state is a localized in the sense that there exist a high prob- 
ability to find the two quanta on the same site, but with equal probability at 
any site of the chain. For the chain with a point impurity, we plot in Fig. [3] 
the coefficients of some of the components of the ground state wave function 
for various values of 7i m - The left hand figure shows the coefficients of the 
[20 . . .], [020 ...],... components. As y lm increases, these coefficients start from 
an initial spatially uniform distribution, but then localize around the site of the 
impurity, in this case at £ = 10. At the largest value of 7; m shown, over 60% of 
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Figure 3: QDNLS model with point impurity, 7 = 4, N = 2, / = 19. 



the wave function is in the state [0 . . . 020 . . .], with the two bosons at site 10. 

The right hand figure shows the corresponding coefficients of the [110 ...], [0110 ...],.. . 
components. Again some localization is found as 7i m increases, but this effect 
is much weaker for these components. Other components are even smaller. 

In Fig. 2] we plot the size of the first three components of the type [2] , [1 1] , [101] 
respectively, centred around site 10, as a function of 7i m — 7. The localization 
increases very rapidly with the magnitude of the impurity. Note that there is 
no Anderson-like effect in this case as the harmonic terms are homogeneous. 



3 The twisted chain 

For the twisted chain as shown in Fig. [5^, the only extra parameter is a m j, the 
strength of the long range coupling between the two spatially adjacent sites m 
and I. As an example we consider the case a m .e = 1, / = 19, m = 5, and £ = 15. 
Fig - El shows two sets of components of the ground state wave function, plotted 
as a function of 7. The components here are the same as those plotted in Fig. 
with a crossover point at m = 5, £ = 15. A breather localized at the crossover 
point will show an enhanced coefficient corresponding to localization at the 
two points of the chain which come together. With 7 = 0, we see some small 
localization effect at sites m and £ for the [2] coefficients at the crossover points, 
corresponding to a harmonic Anderson- like effect. However for nonzero 7 we see 
that this effect is strongly enhanced. The coefficient of the [11] component show 
a weaker localization as before. Similar results are obtained for other values of 
otmi 1 with the strength of the localization depending on the size of a m ,i- These 
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Figure 4: Some components of the wave function corresponding to the ground 
state, N = 2, / = 19, 7 = 4. 



trapped quantum states simulate the trapping of the classical mobile breather 
studied in 0. 

4 The Bent Chain 

For the bent QDNLS chain Fig. we follow the classical DNLS treatment \12\ . 
We consider the Hamiltonian jljona finite lattice with the additional term (J2J) . 
The parameter a is related to the wedge angle 9 through at = k{l — cos#) _1 , 
and we take the site of the vertex to be mo = \{f + 1). 

Fig. shows the coefficients of the [2], [11], [101] components of the wave 
function for various values of 9. Fig. shows some components of the ground 
state wave function corresponding to the neighbors of the vertex, for both 7 = 
and 7 = 4. If the angle 9 is close to 7r, the behavior is similar to the straight 
chain. The ground state is weakly localized around the center (vertex) of the 
chain. As 9 decreases, the localization around the vertex increases, and when 
this angle is small enough, the largest components of the wave function in the 
ground state consists of states localized around the vertex and the two connected 
neighboring sites. In the limit 9 — > 0, the lattice becomes a T-junction, a 
model of interest in its own right. It is interesting that the localization in the 
anharmonic model exhibits a maximum at 9 k, 0.5, whereas in the harmonic 
case the maximum is at 9 = 0. Also the enhancement due to the anharmonic 
terms goes to zero as 9 — > 0. 

In this paper, we have considered only a long-range interaction between the 
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Figure 5: QDNLS model with twisted chain, 7 = 4.0, N = 2, f = 19 
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Figure 6: QDNLS bent chain model, N = 2, / = 19, 7 = 4. 
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Figure 7: QNLS bent chain, N = 2, / = 19, various 7. 

two vertices of the chain. We have also studied more realistic models with long- 
range interaction between all neighbors in the chain. Qualitatively the same 
localization phenomena is observed: a full description will be given elsewhere. 
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